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ON THE TRINOMIAL UNILATERAL QUADRATIC EQUATION IN 
MATRICES OF THE SECOND ORDER. 


[Quarterly Journal of Mathematics, xx. (1885), pp. 805—312.] 


In the May number [p. 225 above] of the present year of the London 
and Edinburgh Philosophical Magazine (disfigured by numerous errors or 
inaccuracies) I investigated the number of the solutions of an equation in 
quaternions or matrices of the second order, belonging to what I term the 
unilateral class, meaning one in which the coefficients of any actual power 
of the unknown quantity lie on the same side of it; this number for the 
Jerrardian Trinomial form I proved strictly is 27? — i (i being the degree 
of the equation) and with evidence little short of moral certainty 7 — i + 1* 
in the general case where none of the terms are wanting t. 


But it must be well borne in mind that these numbers only apply when 
the coefficients are left general, and that for special relations between them 
some or all of the roots may become either ideal or indeterminate, or some 
the one and some the other. In all cases of equations in matrices one 
principal feature of the investigation is, or should be, to determine the 
equation of condition between the coefficients, in order that the solution 
may lose or retain its normal form; if we wish to avoid being compelled 
to enter upon a complicated consideration of exceptions piled upon ex- 
ceptions, it is necessary to presuppose a certain criterion function to be 
other than zero; otherwise it is like the opening of Pandora’s box, letting 
loose an almost incalculable train of vexatious inquiries scarcely worth. the 
trouble they give to answer correctly. 


* This article was written and sent to the press many months ago. I have since shown that 
the number of roots of a general unilateral equation of degree i in matrices of the order w is the 
number of combinations of iw things taken w and w together, and consequently for the case of 
quaternions is 2i?—i for the general and not merely for the Jerrardian form. See [above, 
pp. 197, 233. Also] Nature, Nov. 13, 1884. 

+ Tmade the assumption that the required number is an analytical function of w. 
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Take as an instance the subject of monothetic equations. I have defined 
a monothetic equation to be one in which all the coefficients are functions of 
a single matrix, which may be called the base. In such an equation of the 
degree i and of the order œ in the matrices, we may suppose the unknown 
quantity to be a function of the base, and then the general formula for 
expressing a function of a matrix as a rational and integral function of the 
matrix with the aid of its latent roots, shows that 7 and no more of such 
roots exist. But this in no manner precludes the possibility of the existence 
of other roots which are not functions of the base. Thus, for example, in the very 
simple case of the equation æ? + px = 0, where æ and p are quaternions or 
matrices of the second order, I have shown in the Comptes Rendus [pp, 174, 
179 above] that besides the four determinate ones, all of which (0 included) 
may be regarded as functions of p, there are two other indeterminate ones, 
each one containing an arbitrary constant, and neither of them (to use 
quaternion language) coplanar with the base. Here there is a sort of rever- 
sion to the normal case of 3 pairs of roots to an unilateral quadratic, with 
the modification of two of them having become indeterminate. It becomes 
then of importance to fix accurately the condition of this normal state of 
things ceasing to exist. Being intent on the Denumeration theory of the 
roots in the general quadratic, I did not in the paper cited do this explicitly 
for the unilateral quadratic, although I gave there my own form of solution. 
Moreover, there are other features of much interest belonging to the 
question, which, for the same reason, I omitted to notice. These omissions 
and shortcomings it is the object of this present article to supply. 


Starting with the form a? — 2ps + q = 0, and for convenience of com- 
parison with Hamilton’s formule treating p, q indifferently as matrices or as 
quaternions, and forming the equation 4? — 2Be + D=0, where B, D are 
scalars to be determined, so that B= Sz and D = Tz’, we shall have 

2x =(p— B)*(q—D). 
If now we understand by b, c, d, e, f 
Sp, Sq, Tp, S( Vp Vq), Tq respectively, 
by means of the general formula 
Tr. (myx? — 28 (Var Vy) (rx) + Tx? = 0%, 

[remembering that 

T (p—By=d? — 2bB + P, 

T (q — DY = f? —2cD + D*, 

* This formula, which I have not met with in Treatises on Quaternions, is a particular case 

only of the general Theorem in Matrices, that if 


AK + BNO wh... + Dp” 
is the determinant to (AL+M), where L and M are two matrices of the order w and \ and p 
two ordinary quantities, then 
A(L71M)® - B(L-1M)” 1... +(-)°L=0. 
8. IV. 18 
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and S {V (p — B) V (9 — D)} = e — bD — cB + BD), 
we shall obtain [see p. 188 above] 
4 (d — 2bB + B?) 2 — 4 (e — bD — cB + BD) « + (f — 2D + D*)=0. 
Hence, writing B-b=u, D-c¢c=v, 
d- =a, e-be=8, f-C=y, 
and comparing with each other the two quadratic equations in #, we may 


write 
W+a=r, uw+hB=2rA(ut+b), V+y=4rA(v +0). 


Eliminating v from the latter two equations there results 
— (2rAu + 2brA — BF + 4A (2AUu + 2b — B) u — (y — 4cr) u? = 0, 
and finally writing A — a for u’, we obtain 
(442 + 4cr — y) (A— a) — (20 — By? = 0. 

There are thus 3 pairs of roots, for to each of the three values of à corre- 
spond two values of u, namely. 
+(~1—d +0), 
and to each value of à and u one value of v, namely 

Or + 2br + be Tai 
u 

We have also æ=+4 [(p—b— u) (q — D), 
consequently, since p? — 2bp + d = 0, 
p—b+u)(y—c—v)_ (p—b+u)(q-c-2) 

2 (b — d — vw’) 2A j 

Thus then we see that æ can only cease to have 6 determinate values 
when à =0, and consequently the Criterion of Normality is the last term in 
the equation to ^. 


HAR 


This equation, written out at length, is 
AN + 4 (c — b? — a)r + (— 4ca + 408 — y) + ay — B? =0, 
that is, 4A? + 4(c —d) A? + (— 4ed + 4be— f + e)a + (d— b) (f—) —(e— bey 
Hence the Criterion in question is (d — b) (f — œe) -— (e — be? or 
df— ed —bf— & + 2bce, which is the discriminant to the quadratic form 
A? + Drv + 2cuv + du? + 2euy + fr’; 
this, as I have elsewhere shown, is the Criterion of the matrices p, g* being in 
involution t, that is, of a linear equation existing between the matrices 1, p, q, 
pq; or if p, q are regarded as quaternions, it is the condition of the square of 
* When p, q are regarded as matrices, then 
p?-2bp+d=0, q°-2cq+f=0, 3 (pqt+gp)—bg-cp+e=0, 
where A? + QAv + 2cuv + du? + 2euv + fv? 


is the determinant to A + up +vq. 
[t+ Above p. 116.] 
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the sine of the angle between the vectors of p and q vanishing; a condition 
which of course does not imply the coincidence of the vectors unless accom- 
panied by the futile limitation of such vectors being real. 


It admits of easy demonstration by virtue of the foregoing that in the 
case of the more general equation 
pe+qr+r=0, 
the Criterion of Normality will be the discriminant of the ternary quadratic, 
which is the determinant of 
putqtrrw; 
this seems to me a very remarkable and noteworthy theorem. When this 
Criterion does not vanish, the quadratic equation above written must have 
3 pairs of determinate roots. 


Why they go in pairs and can be found by solving only a cubic instead of 
a sextic is best seen à priori by reverting to the original form æ? — 2ps +q=0. 


It follows from the nature of the process for finding B and D that they 
will be the same for that equation as for the equation 4? — 2yp +q=0. 


But on writing æ + y= 2p these two equations pass into one another. 


Hence each value of B, say B,, will be associated with another value, say 
B’, where B,+ B’ =2b*, that is to say, if u, namely B —b is one value of u, 
then b— B, that is, — u will be another value of u, so that the equation in «u? 
ought to be (as it has been shown to be) a cubic. 


It might for a moment be supposed that A= a= d — b? would lead to a 


i 5 2br + bc — 
breach of normality on account of the equation v — 2) = DAE oes , where 


wv=0. 
This, however, is not the case. For the equation 
v+y=4rA(u +c) 
becomes, when A= a, 
v—4(d—b?)v+ f- e- 4cd + 4b’ = 0, 
so that v remains finite; consequently 2b + be —e, that is, 2bd — 2b* + be — e, 


must vanish when A = d — b’, and v — 2) assumes the form 2. Obviously then 
in this case, to the one value w=0 will be associated the two values of v, say 
v and v,, given by the above quadratic, and to A = a will still correspond two 
values of (u, v), namely (0, v), (0, v); where, ideally speaking, the two zeros 
may be regarded as the same infinitesimal affected with opposite signs. 


* In quaternion phrase, if £ +y=2p, Sx+Sy=2Sp. 

It should be observed, in order to understand what follows in the text, that b- B,=B’-b, 
and that the values of B must obviously be the same in the equation 2?-2pxr-+q=0 as in the 
equation x? — Qap + q=0. 

18—2 
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The equation in à may be made to undergo a useful linear transformation. 
Let A= u + a, so that w= u. 
Then 
p 4p? + (8a + 4c) u + 4a? + 4ca — y} — (2bu + 2ba — BYP = 0, 
that is 4u? + {4 (c + 2d) — 12°} u? + {(e + 2d)? — 8 (c + 2d) D° 
+ 12b*+ 4be— f } w— {b (c + 2d) — 2b? — e}? = 0, 
where it is noticeable that the number of parameters is reduced from 5 to 4, 


cand d only appearing together in the linear combination c+ 2d. This is 
tantamount to the form obtained by Hamilton. 


Let us make another linear transformation suggested by the preceding 

remark. Write c+ 2d=g, and »—lb?=y=d-—d, the equation becomes 
Any? + Agr? + (9g? + 4be— f) y + 2beg — Pf -—e =0. 

But obviously, notwithstanding this reduction of the parameters, ) itself 
is the most natural quantity to employ as the base of the solution, or, so to 
say, as the independent variable, and this admits of being determined by an 
equation of extraordinary simplicity. 

For, let J be the discriminant of 

det. (A+ pp + vq) =L = df + 2bce — ed —bf — e. 

Then it will be seen by actual inspection that the equation found for ^ 

takes the following form 
A-d T= 0, 
that is 


T+(2 


J I+ 


og alee ale ae Wee at ae ee 


d 2) a (24 d 


1 (24 d 7.x =0, 


(the terms in the exponential function subsequent to the fourth term adding 
nothing to the value of the series). \ 


If in the equation 2° — 2pæ + q =0, p and q be regarded as quaternions, 
then A= Sæ + Ip — (Sp), c= Sq, d= Ip, and I=4(pq— gp), which is a 
scalar quantity, and is to be regarded as an explicit function of Sp, Sq; Tp, 
S(VpVq), Tq’; it is in fact the discriminant of the form 


X? + 2SpX Y +28qXZ4+ Tp? VY? +28 (VpVq) YZ + TEZ, 
an identity unknown I believe to the geometrical quaternionists. 
[As an example of it, let p = ù, q =j, then 
Sp=0, Sq=0, S(VpVq)=0, Tp=—1, Tq=-—1, 
4 (pq — gp)? = 1 = the discriminant of X? — Y? — Z>] 


With these definitions eò (?%- ĉa) T becomes identically zero. 
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The equation a? — 2pæ + q = 0 having six roots it is natural to inquire as 
to the value of their sum. This may be readily found as follows. We have 
found 

__(p-b+u (g-e) 


2N 
Also, if æ+ a= 2p, 
a? — 2x'p +q =0, 
and obviously Sas Sa. 
Hence lei, saat Seem’ 
2A 
and a e 1 E E 
2A 
Therefore Lae = 6p — > ge (pq — gp) 
1 
=Go +3733, — 
6p — 32> = 


(98 — 
mop +32 Oh 84) I 


= 6 {p + (28, — èa) I}, 
where the sign of ZŻ must be so taken that it shall be equal to $ (pq — qp). 
So again La? = 2pLa — ôq 
| = 12p? — 6q + 12 (28, — 84) IŻ p. 

Thus the mean value of each root is e in excess, and that of each square 
root ep in excess, of what these means would be if p and q were nominal 
quantities, e denoting (28, — 84) Ip. Of course Sv’ may be found by the 
formula of derivation 

) Lat! = 2p Lat — 92a‘, 

In conclusion it may be observed in regard to the equation 4? — 2px + q = 0, 
(since in writing æ + æ= 2p, we have a, —2a,p+q=0) it follows that (what- 
ever be the order of the quantities p and q) the roots of either equation must 
be associated in pairs; because, if the identical equation to p is p*—wbp*+... 
and to æ is a —wBa*+..., the equation for finding B must be of the form 
T(B-by=0. 

P.S.—Since the above was sent to press I have discovered the general 
solution of the unilateral equation of any degree in matrices of any order ; 
see the Comptes Rendus of the Institute for Oct. 20, 1884 [pp. 197, 233 
above], and Nature for Nov. 13, 1884*. 


[* This paper contains the Theorem “‘ Every latent root of every root of a given unilateral 
function in matrices of any order, is an algebraical root of the determinant of that function 
taken as if the unknown were an ordinary quantity, and conversely every algebraical root of the 
determinant so taken is a latent root of one of the roots of the given function.’’] 
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